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I. INTRODUCTION 


Recently there has been remarkable progress in the study of exclusive charmless Bu,d 
decays. Experimentally, many two-body non-leptonic charmless decays have been ob¬ 
served by CLEO and R-factories at KEK and SLAG (see Refs. HI], p|, p|, [5|, |6|, |7|, ^ p|, p!0[| ), 

and more B decay channels will be measured with great precision soon. With the ac¬ 
cumulation of data, the Standard Model can be tested in more detail. Theoretically, 
several attractive methods have been proposed to study the nonfactorizable effects in 
hadronic matrix elements from Erst principles, such as QCD factorization (QCDF) 0, 
perturbative QCD method (PQCD) 0 Q, and so on. Intensive investigations on 


hadronic charmless two-body B^^d decays have been studied in detail, for example, in Refs. 

0, n, 0, m n, mmm - 

The potential Bg decay modes permit us to overconstrain the unitarity CKM matrix. 
This makes the search for CP violation in the Bg decays highly interesting. The problem is 
that Bg mesons oscillate at a high frequency, and nonleptonic Bg decays have still remained 
elusive from observation. Today only some weak upper limits on branching ratios of several 


charmless hadronic decays are available, mostly from LEP and SLD experiments |^, such 
as Bg —>■ TT+TT”, vr'^TT^, T]7r^, rji], K~^K~, ■■■. Unlike Bu,d mesons, the heavier Bg 

mesons cannot be studied at the R-factories operating at the T(4S') resonance. However 
it is believed that in the future at hadron colliders, such as CDF, DO, HERA-B, BTeV, 
and LHCb, the signs of CP violation in Bg system can be observed with high accuracy in 
addition to studies of certain B^ d modes . 


The early theoretical studies of two-body charmless nonleptonic decays of Bg meson can 
be found in Refs. 0, 0, 0, 0, 0]. The investigation on the exclusive charmless Bg 
decays into hnal states containing meson was given within the generalized factoriza¬ 
tion framework |0]. Chen, Cheng, and Tseng calculated carefully the branching ratios for 
charmless decays Bg —>■ PP, PV, VV (here P and V denote pseudoscalar and vector mesons 
respectively) [0|. And new physics effects in Bg decays was considered in |0. It is found 
that the electroweak penguin contributions can be large for some decays modes 0, 0 , and 
that branching ratios for Bg rjrj' and several other decay modes can be as large as 10“® 
0 which is measurable at future experiments. 

Few years ago, Beneke, Buchalla, Neubert, and Sachrajda gave a QCDF formula to 
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compute the hadronic matrix elements {MiM 2 \Oi\B) in the heavy quark limit, so that the 
hadronic uncertainties enter only at the level of power corrections of 1/mb- This basic 
formula is presumed to be valid for B decays into two light hnal states 0,0. We made 
a comprehensive analysis on exclusive hadronic B^^d decay using the QCDF approach, and 
calculated the branching ratios and CP asymmetries for decays B^^d PP and PV p2 


We hnd that with appropriate parameters, most of our predictions are in agreement with 
the present experimental data. In this paper, we would like to apply the QCDF approach 
to the case of Bg mesons. 

This paper is organized as follow: In section |1|, we discuss the theoretical framework 
and dehne the relevant matrix elements for Bg PP, PV decays. In section B we 
list the theoretical input parameters used in our analysis. Section EE and section 0 are 
devoted to the numerical results and some remarks of CP averaged branching ratios and 
CP-violating asymmetries, respectively. In the mean time, the theoretical uncertainties 
due to the variation of inputs are investigated. In section 0, we give the values of the 
penguin-to-tree ratio P-k-k/T-k-k and a constraint on weak angle 7 . Finally, we conclude with 
a summary in section ivi. 

II. THEORETICAL FRAMEWORK FOR B DECAYS 

A. The effective Hamiltonian 


Using the operator product expansion and renormalization group equation, the low energy 
effective Hamiltonian relevant to nonleptonic B decays can be written as [0^] : 

V‘^q=u,c k=3 

+ + C^gQ^g^ + Y1.C.., (1) 

where Vq = Vqb V*^ (for b^d transition) or Vg = Vqt V*g (for b^s transition) are CKM 
factors. are Wilson coefficients which have been reliably evaluated to the next-to- 

leading logarithmic order. Their numerical values in the naive dimensional regularization 
scheme at three different scales are listed in Table |. The effective operators, Qi, can be 
expressed explicitly as follows: 


Ql — {Uaba)v-A{<lpUi3)v-A, 


Ql — {Caba)v-A{(ll3C^)v-A, 


( 2 ) 
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Q 2 = iUahp)v-A{,q(iUa)v-A, 

Q 2 = {Cahp)v-A{qpCa)v-A, 

(3) 

Q 3 = {Qaba)v-A^{d'f3q'i3)v-A, 

9 ' 

Qa = {dl3ba)v-A^{qa(l'f3)v-A, 
q' 

(4) 

Q 5 = {Qaba)v-A^{q'f3q'i3)v+A, 
q' 

Qe = {dl3ba)v-A^{qa(l'f3)v+A, 
q' 

(5) 

3 

Qy = 2^Qaba)v-A^eqi{q'pq'p)v+A, 
q' 

3 

Qs = 2(^pba)v-A'^eqf{q^qp)v+A, 

q' 

( 6 ) 

3 

Qo = ^idaba)v-Aj2^q'(^Wv-A, 
q' 

3 

Qio = 2 ^Qpba)v-A'^eq/{q'^q'p)v-A, 

q' 

(7) 

g 

Q 77 

Qsg - ^^^mbqa(J^’'{l + 

( 8 ) 


where q' denotes all the active quarks at scale /i = 0{mb), i.e. g' = u, d, s, c, b. 


B. Hadronic matrix elements within the QCDF framework 


To get the decay amplitudes, the most difficult theoretical work is to compute the hadronic 
matrix elements of the effective operators, i.e. {MiM 2 \Oi\B). Phenomenologically, these 
hadronic matrix elements are usually parameterized into the product of the decay constants 
and the transition form factors based on the naive factorization scheme (NF) |^. However, 
one main defect of the rough NF approach is that hadronic matrix elements cannot make 
compensation for the renormalization scheme- and scale- dependence of Wilson coefficients, 
in this sense NF’s results are unphysical. This indicates that “nonfactorizable” contributions 
from high order corrections to the hadronic matrix elements must be taken into account. 

The QCDF approach is one of novel methods to evaluate these hadronic matrix elements 
relevant to B decays systematically. In the heavy quark limit mb 3> Aqcd, up to power 


corrections of order of AqcD/fnb, the basic QCDF formula is fn 


+ [ [ dx [ dy 

JO Jo Jo 

= (M1M2I 1 -|- + OjAgcp/mb) 


(9) 


where denote hard-scattering kernels. At leading order, T/ = 1, T/^ = 0, the QCDF 
formula (P) shows that there is no long-distance interaction between M 2 meson and {BMi) 
system, and reproduces the NF’s results. Neglecting the power corrections of O^AQco/^b), 
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T/’^^ are hard gluon exchange dominant, and therefore calculable order by order with per¬ 
turbative theory. Nonperturbative effects are either suppressed by l/rub or parameterized 
in terms of mesons decay constants, form factors ^ and meson light-cone distribution 

amplitudes $^(0) ^m{.x). The factorized matrix elements {MiM 2 \Ji®J 2 \B)p is the same 
as the dehnition of the BSW approximation |^. Through the QCDF formula, the hadronic 
matrix elements can be separated into short-distance part and long-distance part, and the 
“residual” renormalization scheme- and scale- dependence of hadronic matrix elements could 
be extracted to cancel those of the corresponding Wilson coefficients, so that physical re¬ 
sults at least at the order of as level are renormalization scheme and scale independent 
3^. Through the QCDF formula, “nonfactorizable” effects can be evaluated, and partial 


information about the strong phases can be obtained. 

It is important to note that some power suppression might fail in some cases because the 
b quark mass is not asymptotically large. For example, power correction proportional to 
2m\j/{mhmq) with q = u, d, s, which is formally power suppressed, is now chirally enhanced 
and numerically important to penguin dominated B rare decays. Therefore it is necessary 
to include at least the chirally enhanced corrections consistently for phenomenological ap¬ 
plication of QCDF in B decays. However, the twist-3 corrections to hard scattering kernels 
cannot provide sufficient endpoint suppression, so there appears infrared logarithmic 
divergence, Jdx/x ~ \n{mb/A qqd)- In PQCD method, this singularity can be smoothed 
out by introducing the partonic intrinsic transverse momentum and the mechanism of Su- 
dakov suppression. It should be interesting to investigate the possibility of incorporating 
the Sudakov form factor into the QCDF approach [^. But it is really uneasy because 
the Sudakov suppression is one of the key ideas of the PQCD method while PQCD and 
QCDF have different power expansions which lead to completely different understanding on 
B decays. Therefore, to take the chirally enhanced corrections into account, in this paper 
we adopt phenomenological treatment for the divergent integral Jdx/x [H- 


For weak annihilation contributions, they are believed to be very small with the naive 


factorization assumption (see, for example Ref. |^|). Within the QCDF approach, the weak 
annihilation amplitudes are also formerly suppressed by ~ Agco/'^b 

Q. But as emphasized in the PQCD method |^, |^, annihilation contributions with 
QCD corrections could give potentially large strong phases, hence large CP violation could 
be expected. In addition, the phenomenological investigations on B decays 0 0 0 
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within the QCDF framework also suggest that their effects could be sizable when large model 
uncertainties are considered. So annihilation contributions cannot be simply neglected. In 
the previous work ^ it has been shown that annihilation contributions exhibit 


endpoint singularities even with leading twist distribution amplitudes for the hnal states, 
and these infrared divergence must be parameterized, so extra theoretical uncertainties and 
model dependence are introduced. In spite of these problems, it is still interesting to estimate 
the weak annihilation effects in this paper. 

With the QCDF approach, matrix elements for two-body Bg decays can be written as. 


Gf 


10 


V ^ q=u,c i=l 


G 


A^iBg^M,M 2 ) oc ^ E Y.fBjMjM,V,h, 

V ^ Q=U,C i 


( 10 ) 
( 11 ) 

( 12 ) 

q=u,c i 

where A°‘ arises from weak annihilation contributions, /m are decay constants for 5^, 
and M mesons, respectively. (M can be either pseudoscalar meson P or vector meson V) 
The explicit expressions of decay amplitudes A-^'°' for Bg —> PP, PV are listed in Appendix 
01 ,00. Summary of the dynamical quantities ai, bi is given in Refs, 


III. INPUT PARAMETERS 


Theoretical expressions for decay amplitudes using the QCDF approach are complicated 
and depend on many input parameters including the SM parameters (such as CKM matrix 
elements, quark masses), Wilson coefficients and the renormalization scale /i, and some soft 
and nonperturbative hadronic quantities (such as meson decay constants, form factors, and 
meson light cone distribution amplitudes), and so on. If quantitative predictions are to be 
made, the values for various parameters employed in this paper must be specihed. It has been 
shown that the renormalization scale dependence has been greatly reduced compared to the 
NF’s coefficients obtained at leading order |2^, |^, and the residual scale dependence 
should be further reduced when the higher order radiative corrections are considered. In 
calculations we use fi = mb. The rest parameters are discussed below. 
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A. The CKM matrix elements 


VcKM — 


+ 0 (\‘‘ 


(13) 


One widely used approximate form of the CKM matrix is the Wolfenstein parameteriza¬ 
tion [^1 which emphasizes the hierarchy among its elements and expresses them in terms 
of powers of A = |14s|, 

^ 1 — A ^/2 A A\^{p — ir])^ 

-A 1 - AV 2 AA^ 

AX^ {1 — p — ip) —AA^ 1 

The values of four Wolfenstein parameters (A, A, p, and p) are given by several analysis 
methods from the best knowledge of the experimental and theoretical inputs (for example, 
see Table 0 . Within one standard deviation, the results from different approaches |^, 
n| , |i2[| are virtually consistent with each other. In this paper, we shall take A = 0.2236 ± 
0.0031, A = 0.824 ± 0.046, p= 0.22 ± 0.10, p= 0.35 ± 0.05, and 7 = (59±13)° [g3 . 




However, it is not the right time to draw dehnite conclusions on the parameters p, p 
and 7 . Some interesting hints seem to favour 7^90°, which is in conflict with the data in 
Table ||. For example, it is assumed that it is possible to derive constraint on the angle 7 
from a global analysis B decays. Bargiotti, et al. obtain the bound I 7 — 90°| > 21° at 95% 
C.L. from B —y Kti decay rates and CP asymmetries H- Combining the QCDF approach 
with a global CKM matrix analysis — Rht scenario advocated in |^, Beneke and Neubert 
make a £t of the Wolfenstein parameters {p,p) to six B —> tttt, Kn decays with updated 
measurements, and their fit tend to favor 7 > 90° 1^^. In analogy with works |Q, we also 
make a fit of the CKM matrix parameters to B ^ PP, PV decays, and the preliminary 
results are A = 0.22, A = 0.82, p = 0.086, p = 0.39, and 7 = 78.8° |^. For comparison, we 
shall take the results in as the CKM matrix inputs. 


B. Quark masses 

There are two different classes of quark masses. One type is pole mass for constituent 
quark, which appears in the penguin loop corrections with the functions Gm(S(j) and 
where Sq = my ml. The definitions of Gu^Sq) and GM{,Sq) can be found in |^. In this 
paper, we take 

= md = ms = t), mc = 1.47 GeV, mb = 4.66 GeV. (14) 
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The other is current quark mass which appears in the equations of motions, and is renor¬ 


malization scale dependent. Their values are p3 

1 


m„(2GeV)+mrf(2GeV)J = (4.2±1.0) MeV, 
m,(2GeV) = (105±25) MeV, 
mb{mb) = (4.26±0.15±0.15) GeV. 


(16) 

(16) 

(17) 


Here we would like to use their central values for discussion. And using the renormalization 
group equation, their corresponding values at the scale of /i = 0{mb) can be obtained. 
Because the current masses of light quarks are determined with large uncertainties, for 
illustration, we take approximation 


2 /ip 


TT 71 

= r = r' 
rUb ^ ^ 


(') 


fu 

/„(') 


2ni\ 


= r^ = ■ __ 

^ mb{ms + rriq 


(18) 


C. Nonperturbative hadronic quantities 


Nonperturbative hadronic quantities, such as meson decay constants, form factors, and 
meson light cone distribution amplitudes, appear as inputs in the QGDF formula (|^). In 
principle, information about decay constants and form factors can be determined form ex¬ 
periments and/or theoretical estimations. Now we specify these parameters. In this paper, 
we assume ideal mixing between uj and cp, i.e. u = {uu + dd)/\/2 and (p = ss. As to r] and r]', 
we take the convention in |^, 46], using two-mixing-angle formula for the decay constants, 
but without the charm quark content in r] and rj'. 


(0|g7/i75g|h^'^(p)) = {q = u, d, s) 

(0|h75M|7/W) /“(O 


m%) 


r 


(0 


(0|s75s|?7W) 

/“ = 


(0|s75s|r/W) = 


And for Bg 

T^Bs-T] _ 

Bq — 


(19) 

( 20 ) 

( 21 ) 

( 22 ) 


transition form factors, we take 

^cose + ^sin0)F(f“"-^ = ( - ^sine + -LcosO) 

VE Vs ^ ^ ^ Vg Vs ^ 

The values of these parameters are collected in Table B 


(23) 


















In this paper, we consider the contributions from chirally enhanced twist-3 light cone 
distribution amplitudes of a light pseudoscalar meson. As to vector mesons, only the longi¬ 
tudinally polarized twist-2 terms are taken into account, and the effects from transversely 
polarized and higher twist parts are neglected because they are power suppressed. In calcu¬ 
lation, we shall take their asymptotic forms, as displayed in |^, i.e. for a light pseudoscalar 
meson, we have [^, 


{P{k)\q{z 2 )q{zi)\Q) 

ifp 


twist-2 asymptotic forms: <hp(a;) = 6xx, 
twist-3 asymptotic forms: ^^p{x) = 1, ^p{x) = 6xx. 

where fp is a decay constant; z = Z 2 — Zi, x = 1 — x. 

For a longitudinally polarized vector meson, we have |^, 

nl 


(24) 

(25) 


(0|g(0)7/,g(z)|K(/c, A)) = k^ 


e^-z 
‘ k-z 


dx e 


-ixk-z 


<l>^(a:). 


twist-2 asymptotic forms: $|/(x) = 6xx. 

where e is a polarization vector, and ey = k/mv- 

For the wave function of Bg meson, we take 00 0]: 

mie u:l9^ 


(26) 

(27) 




2ojI 2 


(28) 


where Np is the normalization constant. ^^(O is peaked around ^ ~ 0.1 with ujb = 0.4 GeV 
and 6 = 0. 

As to the divergent endpoint integral / dx/x^ in analogy with the treatment in works 
we parameterize it as 


A = 


"i dx 


= {l + Qp^) In 


mb 


^>< 1 , 


j Q X Ah 

In numerical calculation, we take their default values as 


O°<0<36O" 


(29) 


decay modes 

Qh 

(t>H 

Qa is 

(t>A is 

Ah 20 

Bg^PP 

Bg^PV 

0 

0 

0 

0 

0.5 

1.0 

10 ° 

330° 

0.5 GeV 

0.5 GeV 


where {qh, 4)h) and {qa, 4>a) are related to the contributions from hard spectator scattering 
and weak annihilations, respectively. 
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IV. BRANCHING RATIOS 


The branching ratios for charmless Bg ■ 
written as: 

BR{Bg^M^M2) = 

where 

\P\ = 


PP, PV decays in Bg meson rest frame can be 
rBs \p\ 


Stt ml 


■jA(Bg^M^M2)j^. 


Bs 


“ ("""i + ™M2)^1|™I. - (niM, - 


2mg^ 


(30) 


(31) 


The lifetime and mass for Bg meson are tb^ = 1.461 ps, and ms^ = 5369.6 MeV [^. And 
since the QCDF approach works in the heavy qnark limit, we take the masses of light mesons 
as zero in the compntation of phase space, then \p\ = mBj2. 

The nnmerical resnlts of CP averaged branching ratios for Bg decays are listed in Table 
rV| and Table |V], which are calcnlated at the scale of /i = mb with two sets of CKM matrix 
parameters. The data in Si? colnmns are calcnlated within the NF framework and BP oc 
the data in BP^ and BP^^°‘ colnmns are estimated with the QCDF approach, and 
BR^ (X. BR^^°‘ cx \A^ + In the following there are some remarks 


• Only several interesting decay modes, snch as Bg 

have large branching ratios, which might be observed potentially in the near fntnre. 
Branching ratios of other decay modes are small, not exceeding 1x10“®. Especially 
for decays Bg 'np^'\ 7 r 0 , pri^'\ whose tree contribntions are snppressed by 

both CKM factor and color, and pengnin contribntions are electroweak coefficient ag 
dominant, so their CP averaged branching ratios are very small, aronnd 0(10“^). As 
to these pnre weak annihilation decays, snch as Bg —> tttt, vrp, vro;, their branching 
ratios are extremely small, aronnd 0 ( 10 “^). 


For those 6 —*-s transition decay modes, snch as Bg —> K^*^K, their tree contri¬ 

bntions are CKM snppressed, and “nonfactorizable” effects contribnte a large portion 
to pengnin coefficients 04 ^ 6 ; so pengnin contribntions and tree ones are either com¬ 
petitive, or pengnin dominant. So we can see large “nonfactorizable” effects in these 
decays. In addition, the coefficients bi and/or 63 appear in these decay amplitndes, 
and the data in Tables and 0 show that weak annihilation contribntions are size¬ 
able (>50%). For those b^d transition decay modes, snch as Bg —>• 
they are Oi dominant, and the radiative corrections are ag snppressed compared with 
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leading order contributions, so the data in Table and 0 show no large difference 
between the results obtained with the QCDF approach and the NT’s ones. 

• There are hierarchy among some decay modes, such as 

BR(^^^K+K-) > BR0l^K+K*-) > BR(bI^K-K*+), (32) 

BR(^^^K°K^) > BR(bI^K^K*^) > BR(Bl^TfK^^), (33) 

BR(Bl^K+7r-) > BR(Bl^K*+7r-), (34) 


There are two lines of reason for the above relations. One main line of reason is that 
the penguin contributions are important or/and dominant for these decays. Their 
decay amplitudes involve the QCD penguin parameters 04 and in the form of 04 + 
RaQ, where i? > 0, for Bg —> PR decays, and i? = 0 (i? < 0) for Bg PV decays 
with B^P (B^V) transition, respectively, as stated in |^. The other line for the 
second inequality of Eq. (|32|) and Eq. (^ is fK*Fi‘’^^ > ■ The numerical 

data in Tables 0 and 1^ confirm the above relations in general. Here we would like to 
point out that because the weak annihilation parameters b 3 {P,V) = —b 3 {V,P), and 
the combination of &1 + 63 (or 63 + 264 ) is destructive for —> K^K*~ (or and 

constructive for K~K*^ (or K^K*^), the weak annihilation have more effects 

to decays of ^ K-K*+ (K°K*^) than to decays of 5° ^ K+K*- (iF°Z*°), (for 
example, see the data in Table 0) which might be an explanation of why the data 
BR^+^(bI^K-K*+) > BRf+^0l^K+K*-) in column 4 of Table 0 that violate the 
second inequality of Eq. (l3^). 


• It is interesting to note that Bg decays have large branching ratios. In fact, 

their SU{3) counterpart, B^^d Kt]' have been reported to have the largest branching 
ratios among the two-body charmless rare B decays. 


Decay modes 

CLEO 

51 

BABAR 1 

1 

Belle 

BR{Bd^Kgr]')xlO^ 

BR{Bu^K^7]')x10^ 

89ll^±9 

80l^°±7 

42ll?±4 

70+8+5 

55tl^+8 

79l{?+9 


The abnormally large branching functions for B^^d Krj' decays have triggered in¬ 
tense theoretical interests in understanding the special property of meson ri^'\ Several 
mechanisms have been proposed (for example. Refs [^, ^ ^). There are many 
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works devoted to the study of exclusive B decays into two-body final states containing 
such as Refs. It is generally believed that this problem is related to the 

axial anomaly in QCD, but the dynamical details remain unclear. It is now commonly 
believed that maybe there is large coupling between two gluons and rj' which might 
have important contributions for rj' production [Q, It is shown in that the 
contributions of g*g* —>■ g' to the formfactors can give a good explanation for exper¬ 
imental data. Recently, M. Beneke and M. Neubert computed the exclusive 

+ X decays using the QCDF approach |M]. Their novel idea is to consider the 
flavor-singlet amplitudes for producing g^''^ meson from not only a quark-antiquark pair 
but also a pair of gluons. Their analysis including three effects: b^sgg amplitude, 
spectator scattering involving two gluons, and weak annihilation, could qualitatively 
account for the measurements with inputs of specihed values, but with large theoretical 
uncertainties. Their conclusion is that it is the constructive or destructive interference 
of non-singlet penguin amplitudes that is the key factor in explaining the exclusive B 
g^''> -|- X decays. 


From Table |V| and 0, we can see that CP averaged branching ratios for many decay 
modes are stable against the choices of the CKM matrix parameters. Only those 
decays which have large interference between tree contributions and penguin ones, 
such as Bg —> • • •, are sensitive to the choice of the angle 7 . 


Of course, theoretical uncertainties from input parameters (such as the CKM matrix el¬ 
ements, quark masses, form factors, X, and so on) should be taken into account when 
discussing Bg decays, which has been studied in detail in Refs. |^, In Figures and 

we consider the effects of the variation of inputs on the CP averaged branching ratios of 
Bg —>■ K^K~, and tt+tt" decays on the weak phase 7 . In each plot, 

the dashed lines, solid lines, and doted lines give the QCDF’s predictions with default values 
of various input parameters for m<j(2 GeV) = 90 MeV, 105 MeV, and 120 MeV, respectively, 
but keep the ratio of light quark masses hxed, rflq/mg = 4.2/105. For discussion, we also 
vary the form factor by ±10%, i.e. = 0.250 ~ 0.300. It also includes the 

uncertainties from the CKM matrix parameters. From Figures |l] and we can see that 

• There exist sizeable theoretical uncertainties which smear some helpful information on 
the angle 7 and demote the predictive power of the QCDF approach. 
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• The variation of form factors brings very large uncertainties (see the 2nd row of Figure 
1 ^ which in principle could be reduced by the ratios of branching ratios, while the 
uncertainties from Xh which parameterizes the divergent end-point integral in hard 
spectator scattering corrections are very small (see the 3rd row of Figure |^). 

• For tree dominant decay modes (such as the theoretical uncer¬ 

tainties mainly come from the formfactors and CKM matrix inputs. While for those 
decays with large interference between tree and penguin amplitudes, such as Bg 
K~^K~, and penguin dominant decays, such as Bg —> the theoretical uncer¬ 

tainties originate mainly from the variation of light quark masses and parameter Xa 
besides formfactors. 

• Bg ^ 7r^n~ decay is pure annihilation process. Its amplitude is free from transition 
form factors and hard spectator scattering corrections. Hence the dominant theoret¬ 
ical uncertainties would come form the weak annihilation effects, more precisely the 
quantity Xa (see Figure @). Experimentally, it is worth searching for pure annihilation 
processes which may be helpful to learn more about the annihilation mechanism and 
to provide some useful information about hnal states interactions and nonperturbative 
parameters, such as Xa- 


V. CP ASYMMETRIES 


CP-violating asymmetries for Bg decays has been studies in |2^. In this paper, we 
shall evaluate them with the QCDF approach. In principle, the calculation of CP-violating 
asymmetries for Bg are similar with those for B^ decays. Due to flavor-changing interactions, 
B^g and Bg can oscillate into each other with time evolution. The time dependent CP 
asymmetries Acp for Bg decays is dehned as 


r(g;(t)^/) - r(g.°(t)^/) 


Acp{t) = 


(35) 


As discussed previously in Refs. [^, the Bg PP, PV decays can be classihed into 
three cases according to the properties of the hnal states. 


case-I: P° 


but P° 74 /, ;^ ^4 /, for example, B^ K+p , vr K*+, 
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•, the CP-violating asymmetry for these decays is time independent, 


•^cp — 


r(p;^/) - r(pg^/) 
r(:^-./) + r(po^/)' 


(36) 


case- 


.II: B.« ^ (/ = /) 


P°, for example, Bg 


the time- 


integrated CP-violating asymmetry for these decays is 

1 Xg 


Acp = 


a^i = 


1 — I Ac/ 
1 + |A 


CP I 


Oe+e' — 


Acp = 


/T=rO, 


-2Im(Acp) , _ VuV;, 


1 + I Ac/ 


V,*gV,, A{B0{0)^fy 


(37) 


(38) 


where a^' and 0 ^+^/ are direct and mixing-indnced CP-violating asymmetries, respec¬ 
tively. The parameter Xg = AtjibJTbs is considerably large for Bg system, Xg > 19.0 
at 95% C.L. 1^. In our calculation, we shall take the preferred value in the SM Xg ~ 


20 1^. Clearly, CP-violating asymmetry ^cp should be very small because a^i and 
in Eq. (|37| ) are strongly suppressed by 1/x^ and l/xg, respectively. 

case-III: po ^ (/&/) ^ B^, for example, Bg {K^K*^ & {K+K*- & 

K~K*^), (vr+p” & 7r“p+). Analogous to the notations for Bd decays in the time 
dependent decay widths for this case of Bg decays are written as: 

r(p°(t)^/) = ^ 

r(p°(t)^/) = 




3 r_B t 


3-rss* 


(isr+i'f) 

1 -|- ag>cos{AmBA) + ae+e'sin(AmsA) 

(isr+iAt) 

1 — ag'Cos{AmBA) — ae+e'sin(AmB^f) 

(l9l" + |At) 

1 -|- a-^i cos {Am B A) + ae+?'sin(AmB^t) 

(isr+i'f) 

1 — ae/cos{AmBA) — ae+e'sin{AmBA) 


where 


g = ^(P°(0)^/), -g = AiWgiO)^f), 

h = A0i{O)^f), h = AiB^M^f), 


and with q/p = VtgV^l/VygVtb, 


\9?-\h\ 


ipp+|hr 

0-e+e' — 

\h\^-\9\^ 





-2lm[{q/p)x{h/g)] 
l + \h/g\‘^ 

-2lm[{q/p)x{g/h)] 
1 + \g/hy 


(43) 

(44) 


(45) 

(46) 
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Our numerical results of OP-violating asymmetries for Bg PP, PV decays are listed 
in Tables OH. which are calculated with two sets of CKM matrix parameters. In the 
following there are some remarks. 


From Tables [VlTH^, we can see that as expected, due to the large parameter Xg 
suppression, OP-violating asymmetry Acp for these case-II decay modes are indeed 
very small, not exceeding 


• From the QCDF formula Eq. (P), we know that radiative corrections and “nonfac- 
torizable” contributions should be either at the order of ag or power suppressed in 
therefore, the (7P-violating asymmetries Acp for these ai dominant de¬ 
cays, such as Bg —^ are not large because of the small strong phases. 


From the previous Bu^d PP, PV analysis ^|, we know that “nonfactorizable” 
effects contribute a large imaginary part to the coefficients 02 , 4 , 6 - So from Tables 


VI|-pC|, we see that for these decays which have large interference between tree contri¬ 
butions dominated by 02 and QCD penguin ones, and large interference between CKM 
suppressed 02 dominant tree contributions and electroweak penguin ones, there exist 
large direct CP-violating asymmetries, for example, Acp for P° —> K*^riV 

decays, a^/ for Bg —> decays, and a^/ for Bg —»• 7r°0, 

decays. 


• The CP-violating asymmetries for the most decay modes is not keen on the variation 
of the CKM matrix parameters. The large and direct CP-violating asymmetries which 
are sensitive to the angle 7 are only display for some of these decays, such as Acp for 

^ K:*° 7 r°, and a,, for Bg 

• If we assume that mixing phase is negligible, or q/p = 1, it should be convenient 

to determine the weak angle 7 or /3 from measurements of CP-violating asymmetries 
for Bg decays. Unfortunately, the very rapid oscillations are expected due 

to the large mixing parameter Xg in SM, which makes the experimental studies of 
CP violation in Bg meson system difficult, and only bounds on few decays modes 
are given for the moment. In addition, there exist large theoretical uncertainties. 
Within the QCDF approach, the subleading power corrections in l/ruf, are might be as 
important as the radiative corrections numerically because mi, is not inhnitely large. 


15 












So the QCDF approach can only give the order of magnitnde of the CP-violating 
asymmetries, as stated in . 


VI. EXTRACTING WEAK PHASES FROM Bs^KK DECAYS 


It is necessary to test the self-consistency of the CKM description of CP violation throngh 
a variety of processes. One test involves the Pd(t) —vt+tt" decays which are potentially rich 
sonrces of information of both strong and weak phases. Experimentally, BABAR and Belle 
have reported the measnrements of CP-violating asymmetries in Bd{t) —>■ 7r’''7r“ decays. 



BABAR 0 

Belle 0 

Average 

Ptttt 

-0.02±0.34±0.05 

1 91 + 0 . 27 + 0.13 
-‘-•^■‘--0.38-0.16 

0.48±0.35'5 

r' 

^TTTT 

-0.30±0.25±0.04 

-0.94l|];i^±0.09 

-0.54±0.20 


which has triggered high theoretical interest. Theoretically, if we assnme that the pengnin 
amplitndes are zero for Prf(f) —>■ 7r’''7r“ decays, then it is expected to determine the weak angle 
a from the S'tttt = —sin2a. Unfortnnately, this relation is strongly pollnted by pengnin effects. 
The natnre of tree and pengnin amplitndes lead to some indeterminacy in determining the 
weak and strong phases. Hence, it will be an interesting work to investigate the ratio 
of pengnin to tree amplitndes. Using the P-spin symmetry, we can get some additional 
information on the pengnin-to-tree ratio, Pj^Tr/Pn-Tr, from its connterpart Bg —> K^K~ as a 
cross check. To illnstrate, we now describe the expressions for the decay amplitudes of Bg 
K^K~ and B^ 7i~^7r~ as follows: 




A{B —>'7r’^7r ) = iTjrTrle 
A(bI^K^K-) = \T,\e 




Using the SU{3) flavor symmetry, we have 

1 


IT 1 

1 TTTT 1 

in+ii 

171 


IP 1 

1 TTTT 1 

\VchV:,\ 

171 

+ 0 + 7,1 

IP 1 

IP 

= tan^^c-^ 

IT 1 

1 TTTT 1 

17 


AV2 


A 

A 


1 -A 72 ’ 

, 6t — 


+ |P..|e-*''7 

(47) 

- |Pc|e-*'p 

(48) 

St = S^T, 

(49) 

Sp = Sp, 

(50) 

1 

II 

(51) 


® the error includes a scale factor of 1.3 
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(52) 


and with ta 

Pr 


IBsIk 

FBs^K^2 , 


_ -\ycbVcs\{(^4 + Qip + ^x (^6 + Qs) + r^(&3 + 2&4 - ^^ 3 "’ + Ibl^)} 

Tc \VubV*,\{a\ + al + + r^(a^ + a^) + r^( 6 i + 63 + 264 - i&r + ’ 

Compared with Bd —> tt+tt" decay, the contribution of Tc {Pc) for Bg K^K~ decay is 
reduced (enhanced) by tan 6 *c. Of course, the relations of Eq. (|i9|) and Eq. (^OD are affected 


by [/-spin breaking effects, such as factor ; and so on. But SU{3) flavor 


("^Ss ’ 

breaking effects are expected to be very small in Eq. (|^ ) within factorization approach. 
As stated in |^, Q, assuming that mixing phase is negligible, and taking the angle 

(3 as one known input which can be determined from Bd J/'^Ks decays and has been 
tentatively given by BABAR 


|] and Belle the strong phase 6t — 5p and |P 7 r 7 r/T^ 


(or Sf — Sp and \Pc/Tc\) as a function of weak angle 7 or/and a can be determined from 
measurements of 5'^,^ and Ct,.„ (or Skk and Ckk)- And employing the relation Eq. ® 
within the P-spin symmetry, the penguin-to-tree ratio can be overconstrained from Bd —>■ 
7 r''' 7 r“ and Bg K^K~ decays. And using the value of the penguin-to-tree ratios, some 
information on weak phases 7 or/and a can be extracted from the measurements. 


^■K'K 

^ = 
^TTTr 

r* — 

'-^TTTT 

^KK = 

Skk = 
Ckk = 


V*,V,,A{B%0)^7r+7r-) 


= e 


1 + 

IP /T 

|e*^'e+ 

1 + 1 

P IT 1 

^ TTTT / TTTT | 

gjVg-i 7 


— 2 Im(A.n- 7 r) —sin2Q( -|- 2 P^jr/77r7r COS(5'cOs(a — (3) + 

IP IT 

\\m2(3 

1 1 ^nn 

P 1 — 2 P,r7r/77r7r COsPcOs(q( -|-/3) -|- 

P IT 1 

^ TTTT / TTTT | 

2 


1- |A. 
1 + |A„ 


21 PjrTT / Tnn \ sin(5'sin7 


1 + 21 Ptttt /PTTTr I COs 5 'cOS 7 + I Ptttt / Tw 


= 


V,gV,lA{Bl{0)^K+K-) 
V,*gVt,A{B^M^K+K-) 

—2 Im(Ax_ft:) sin27 — 2|Pc/Tc|cosPsin7 


1 - 

|Pc/T, 

jgiVgn 

1 - 1 

Pc/Tc\ 

gjVg-i7 


1 + |AiCKp 1 - 2 |Pc/Tc|cosPcos 7 + |Pc/Tcp’ 
1 - I >^KK p -21 Pc/Tc I sin(5'sin7 


(53) 

(54) 

(55) 

(56) 

(57) 

(58) 


1 + lAxi^P 1 - 2 |Pc/Tc|cosPcos 7 + jPc/Tcp 
Here the numerical values of penguin-to-tree ratios are given in Table Although using 
different derivation and inputs, within the range of one a, the results of penguin-to-tree 
ratio |P,r 7 r/ 7 A- 7 r| that are calculated with Eq. (^) are in agreement with the result of (28.5 
± 5.1 ± 5.7) % [0 which is calculated with Xa = Xp = \n{mj,/A^), and the result of (27.6 
± 6.4) % (including SU{3) breaking effects), but not as large as 0.41 [|^. In addition. 


the value of the strong phase 6' is also consistent with (8.2±3.8)° |^|. The uncertainties of 
penguin-to-tree ratio are shown in Figure which indicates 6' G (—7°, 22°). 
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In addition, if the penguin-to-tree ratios are determined, then it is expected to extract 
some information or give bound on weak angle 7 from the measurements of Bg K^K~, 
decays in the future. Now let’s illustrate this point. The QCD penguin terms for 
these two decays should be equal according to factorization, so the decay amplitude for Bg 
could be written as, 

A(Bl^TfK^)c:^-\Pc\e-^^p, (59) 


There are two approximations in Eq. (|59D, as stated in |^: ( 1 ) The color suppressed terms 
which are proportional to the CKM matrix factor VubV*^ in Eq. ([M]) and Eq. (|CI| ) are 
safely neglected because of \VubV*g/VcbV*g\ ~ 2%. (2) The tiny isospin breaking effects are 
disregarded, and the small difference of electroweak penguin contributions between these 
two decay modes are neglected. The ratio of CP averaged branching ratios is dehned as 

BR(^^^K+K-) + BR{B°^K+K-) 


Rkk — 


BR{bI^K^ K^) + BR{B^^^WK^) 

1 — 21 Pc/^c I cos(5'cos 7 + |Pc/Tc 


\PclT, 


Then we can get a constraint on 7 , 


cos7> 


Pr 




R 


KK 


(60) 


( 61 ) 


In addition, Gronau and Ronsner also gave a bound on 7 from the decays Bg —>■ K^K~, 
without prior knowledge of the penguin-to-tree ratio, sin^y < Rkk 0 • 

From the above discussion, we can see that sufficient measurements of Bg decays in the 
future can resolve ambiguities on the determination of the CKM angles. Here, we estimate 
the bounds on 7 with the QCDF approach. The data in Table ^ indicates that the weak 
angle 7 given in Refs. [^, and the bound from Eq. (^) are consistent with each other, 
which might be tested in future measurements. 


VII. SUMMARY AND CONCLUSION 

In this paper, we calculated the CP averaged branching ratios and UP-violating asym¬ 
metries of two-body charmless hadronic Bg —> PP, PV decays at next-to-leading order in 
ag with the QCDF approach, including “nonfactorizable” corrections, as well as those from 
weak annihilation topologies. We hnd 
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• Only several decays, such as Bg have large branch¬ 

ing ratios, which might be accessible at hadron colliders potentially in the near fu¬ 
ture to allow for detailed phenomenological analysis. The Oi dominant decays Bg 
TT are generally insensitive to the contributions from both “nonfac- 

torizable” effects and weak annihilation, so the numerical predictions of the QCDF 
approach on their CP averaged branching ratios are similar to their NT’s counterparts; 
the main uncertainties come mainly from the CKM matrix parameters and form fac¬ 
tors. And their direct CP-violating asymmetries Acp are only a few percent because 
of their small strong phases. But for the branching ratios of Bg —> K^*^K, pAfji') 
decays, the contributions of the “nonfactorizable” effects and weak annihilation can 
be sizeable. 


Large direct CP-violating asymmetries occur in some decays whose decay amplitudes 
are related to coefficient 02 , because 02 obtains a large imaginary part from “nonfac¬ 
torizable” effects, and some can even reach 80% or so, for example, a{,{K^rf) —85% 


in Table 0. But Acp for case-II CP decays is very small because of large Xg. Of 
course, it is assumed that the contribution of power corrections in 1/mi, are as impor¬ 
tant as radiative ones to strong phase numerically, so the quantitative predictions on 
CP-violating asymmetries of Bg decays with the QCDF approach should not be taken 
too seriously. 


• Although we can overconstrain the penguin-to-tree ratio, |Pjr^/Pn-jr|, and give a bound 
on 7 from Bg decay into charged and neutral kaons, too many input parameters bring 
large theoretical uncertainties, what’s more, experimental studies on Bg decays are 
very limited so far, so it is not the time to extract some useful information on the 
angle 7 from two-body charmless Bg decays. We look forward to future measurements 
and theoretical developments to give some insight into these parameters. 
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APPENDIX A: THE DECAY AMPLITUDES FOR B^^PP 




-aiQ + -Ri ^06 — —ag) |, 


X < 04 


where Ri = 7—;—#-7. 

^ {rns+md)(mf,-md) 


1 3 

+ iyubVud + ^cbVcd) - 04 + -a-io ~ 2 ~ 

— ^ 2 ( 06 — 2 ^ 3 ) |; 


where R2 = 




AUbI^KV^) = [ml, - m;„.) (K^IC, + K,tC) 

x{a4 — 2^10 ^ ^3^6 ~ 

- m|){l43V'i03 

+ (K(feKrd + ^cbVcd) 2(^3 ~ ^s) + <^4 ~ 2*^^° ~ 2^^'^ ~ 

+4'^ (ag - ^Og) (1 - + (03 - 05 + ^07 - ^09) |, 


where R3 = 


(m3+md){mf,-ms) 


, and i?4 ^ = 


{ms+ms){mf,-md)' 


+ ('43K:, + KsK,;) [ - 1(07 - ag)] |. 

+ (^^6^1*5 + KftKi) 2(03 — 05 ^ — 2 (®7 — a- 9 ) 

+ (yubV*g + VcbV*)j O3 - O5 + O4 - -Oio + -O7 

J'qC) ^ ^ 

-^09 + 4'^ (og - ^og) (1 - ^)] |, 
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where 4'^ = . -.. 

+ (yubV*g + VcbV*g^ 2 (as — as) — -(qt — ag) 

fS 11 

+ {yubV*g + VcbV*g"j as — as + 04 — 

J T} L L 

- 2'*9 + -^ 5(06 - ^as') (1 - I 

+ (KfeKl + ^cbVc^ 2 (^as ~ “ 5 ) - 2 (®7 - flg) 

f^/ 1 1 

+ (K(6Krs + VcbVcs) Os — Os + 04 — -Oio + -07 
Jr)' Z Z 

-ios+ fl^a8 - ios)(l - ^)]|. (A6) 

+ (K(6K,d ^ a4 + aio + Re^ag + ag^ |, (A'^) 

-R* = 

+ (l4iC + 14K:) 04 + Oio + -Rr^oe + og^j I, (^8) 

r-i - 

where = 7 —;— S- - 7 . 

' (m„+ms)(mi,—m„) 

APPENDIX B: THE DECAY AMPLITUDES FOR 


AJ{bI^K^K*^) = V2GFfK*Ff^-^mK4e-PKo){v^bV:, + V,bV*)l^aA-^a,o], (Bl) 


AJ(Bl^TfK* 


= V2GFfKA^‘^^*mK^{e-pj^){VubV:^ + KsKl) 

X |a 4 — 2 ®io ~ (®6 ~ 2 *^®) }’ 
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where Qi 


where Q 2 


where 


2m^ 


{ms+ma){mf,+ma)' 




1 3 3 

— 04 + 2 ®io 2 *^^ 


= GF/c.Ff°^'"m^(e-pKo)|K6K>2 

~\~{yubV*fi + VcbV*^ 2(^03 + 05^ + 04 — 2®io 2*^’^ 2*^^ } 


(B3) 


(B4) 


Af{Bl^K^(t)) = V2GFm^{e-pK^){y^iy:^ + V,,V, 


X 




1 


cd 


1 1 

03 + 05 ~ 2 ®'i' ~ 2 *^® 




04 


~ 2 ®io ~ ^ 2(^6 ~ 2 *^®) }’ 


2 m^ „ 
kO 


(ms+md)(mi,+ms) 

.^0rF*0 


= GFUAf^^*"mK*o{e-p^o)i^VubV:,a2 


cd 


+ (l4ilO + Ol+ 

3 1 

+ 2®9 + < 53(06 ~ 2 


1 3 

— 04 + 2 ®io ~ 2 *^^ 


2m? r 


{mu+md){m},+md) ' 

= GF/7r2l^°^'^m^(e-P7ro)|KbK5o2 

+ (K.C + KfeK 


3 3 

~ 2+^ 1" 2^^ 


(B5) 


(B 6 ) 


(B7) 




V2GFmK*o{e-pFiO) 



0 



VubVud(^2 


~^(yubV*^ + VcbV*^ (2(03 — 05) + 04 — 2®io ~ 2*^’^ 
+ ^09 - - ^a 8 )(l - ^)) 

'' 7^(') 

+ (KfeKd + BcbKrf) (03 - 05 + oo? - 0O9) 
+/,,*Ff(K^KTd + (04 - ^Oio) }, 


(B 8 ) 
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where Oj' = 

+ {V^V:, + KrK:) [|(«7 + 09 )]}, (B 9 ) 

+(^ibVus + VcbV^'j 2^a3 + abj + 2 (**7 + ^ 9 'j |i (BIO) 

= V2GFm^(e-p,<„){/“,„<”-*[V7il^*02 

+ (^^6^1*5 + KftKi) ( 2(03 ~ O 5 ) ~ 2^^'^ ~ 

/*/ 1 

+ (KifeKrs + VcbVcs) {^3 — ^5 + O4 — 2®10 

'' 7^(0 

+ 9®7 - 9^9 - <34^(06 - 9 08)(1 - ^~)) 

^ ^ ^ Jr,(') 

(k^K: + KftV;:) (as + as 

1 1 1 \1 /Th N 

+04 — 2®io ~ 2®^ ~ 2*^®/ j ’ (^H) 

= V2GFfpFf^^^mp{e-pK)[VubV:dai 

+ {yubV*d + VcbV*^ (a4 + aio) |, (B12) 

A^(bI^K+K*-) = V2GFfK*F,^"^^mK*{e-pK)[VubV:,a^ 

+ {VubV:, + KoK:) («4 + aio)}, (B13) 

= V 2 GFfKAf^^*mK*{e-PK){VubV:,a, 

F{yubV*g+ VchV*^ 04 + aio — Qs (ae + ag) |, (B14) 

2m? _ 

where Qs = • 

M(Bl^n-K*+) = V2GFfnAf^^^*mK*{e-p^)[v^bV:,a, 

F(VubV*^ + VcbV*^"j 04 + Oio — <56(^6 + ^s) j; (B15) 

where ' 
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APPENDIX C: THE WEAK ANNIHILATION AMPLITUDES FOR B^^PP 

= -i‘^fBji[{y^v:,+v,,v:.)[h<P’ ,K'>)+h(P’ ,K'‘) 

+h(K'‘,P') - ^bT(K°,K°) - ^bT'(K°,K°) - )]|, (Cl) 

= z^fBjKf.{v^v;j + VB,v;,)!^bAw«j{«)-^-br(w\K<>)y (C2) 
A'“)] + /,‘« P) - \b7P‘, l''’)]}. (C3) 

+ (K.»C + Vb,V:^ [2b47, w") + \bT(7, ir«)] |, (C4) 

+ (KtlC, + Ktic) [|6r’(>r". »)">) + lb7(Py)] |, (C5) 

A^P.^pp) = -iPGBfB.{f;,Bf;,B [v^p,bi(P\P) 

+{v„3v;, + Vb,v;,) {2bAP,P) + lb7(P,P))' 

+p,p, (v,„y:. + Vay;.) [b3(P,P) + byf\ P) 

(C6) 

y(ppp) = -i^fBjp:;.{v„y7[bp',’i) + bliri.P 

+(i4i,k:. + [ 2 * 4 ( 1 )'.»)) + 2 * 4 ( 1 ). ()') 

+i*r(>)'.*) + i'>r('* >)')]} 

{VubV;, + I'd, K’.) 1*3(1 )'.'/) + *3(*. t)') 

+* 4 ( 1 )'. I))+* 4 (>), *') - i*r(’A ’)) - ^*r('). ’)') 

-i*r(>/.-))-i*r(-).-)')}. (C 7 ) 
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= -i^fs,Uf^:(y^,v:, + K4VQ){63(i-, A-+) - i6r(i-, A-+)}(C8) 

+ (K.iIC. + K,v;,) [b,{K-,K+) + h(K*, K-) + h(K-, K*) 
-\br(K-,K*) + hfiK*, K-) - ]^bT(K-,K*)\Y (C9) 

^“(B^^ir-,r+) = + {VuX. + l 4 l.V;,) [ 64 (ir+,ir-) 

+ 64 (ir“,x+) + 64“(7r+,ir“) - ^b‘^(TT- (CIO) 

APPENDIX D: THE WEAK ANNIHILATION AMPLITUDES FOR B^^PV 

A‘(P,^K«p«) = ~fB.fKpV^^V:, + V^V:,)!^b,(pO,K«)-^br(p°,K°)}, (Dl) 

X(B^^A'»W) = ^fB.fKf^{V.,tV:, + KtC){63(1*). A“) - i6r(‘^.-ff”)}. (D2) 

= ^fBjKft{VutVZ, + KsK-J |(>s(A», 0) - '^)}. (D3) 

^“(B^^A»A'“) = ^fBjxfK. (KtC + KiK*.) {*3(A‘”, A») + btiK^K'") 

+64(A‘”, A") - itDA’", A") - i()r(A“, A’") - (D4) 

X(B^^7?“a-“) = ^/Bj,f/,f.(v;3lC, + K3V'^)|()s(A°,A-») + ()4(A-“,7?") 

+64(7?“, A’’") - A"") - i()7(A"“,7?°) - i()r(7r“, A"“)|, (D5) 

2 2 2 J 

^•((S^^jr'A-") = ~fBj,fK.{VutV:, + l43V(;)|63(3r“, A-») - A-»)|(D6) 
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+ (t4klC, + V,tVZ) [264(ir”,p") + 264 (p»,i») 

+\bT(^\p’‘) + \t>T(p'‘y)]\ (D7) 

+{v„tv;, + v^v;,) [^br(^\oj) +>r”)]}. (ds) 

= 0, (D9) 

A‘(ff.^n'bK-0) = ^fB.fK- (k.i.v;- 4 + KiV;-,,) [63(>)<''. A'-») 

+ /‘w [*s(A'-“, r,<'>) - r/'l)] |, (DIO) 

X(B^^p''V) = ^fB.Uf‘i„ (K.1.1C. + t4i.K,;) .#■) + i/'>) 

(Dll) 

A‘(gl^rfPp«) = ^^/b.4/J„{kx:,[6i(>)''',p“)+4i(p",i(‘'>); 

+ (K.i.C + KiK*.) P”) + ll>T(p°, P<'>)]}, (D12) 

+ (14i,Ip:4 + 14i,K:) [26i(p'''x) + 264)^, p"’) 

+i6r(p<'>x) + i6r(‘^. >(">)]}. (D13) 

X(B^^A+p-) = ^/b.4/a-(Ki1Ci + KiV'i) {6s(p-, A+) - ^bfip-. A4},(D14) 

X(B^^A+A-) = ^/b,/a-/p-{v4iK:.[''i(A+,A-) 

+ ('7,1,14, + 1417 :) [biiK-, A+) + 64(A+, K-) + btdC-, K+) 
-IbTiK-, A'+) + 6r(A+, A—) - ^bf{K-, A'+)] |, (D15) 
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+ {VubV:, + V^hV*) [h{K-, K*+) + b^{K*+, K-) + h,{K-, K*+) 
-hl^{K-, K*+) + bl^{K*+, K-) - h^^{K-, K*+) 


(D16) 


A^0i^7r+p-) = ^f^J^f^\V^,v:,b^i7r\p-) 


A‘^0i^7r-p+) = 6i(p+,7r 


72 ' 


7r+,p ) + &4(p ,vr+) 


>7r+)] 

[■ 

(D17) 

) 

p+,7r ) 

+ bi{Ti ,p+) 



[. 

(D18) 


l]t)19) 
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079901, (2002). 


TABLE I: Wilson coefficients Ci in the NDR scheme. Input parameters in numerical calculations 
are: as{mz) = 0.117, aemimw) = 1/128, mw = 80.42GeV, mz = 91.188GeV, mt = 178.1 GeV, 
nib = 4-66 GeV. 




= rnb/2 


= mb 


= 2mb 

NLO 

LO 

NLO 

LO 

NLO 

LO 

Cl 

1.130 

1.171 

1.078 

1.111 

1.042 

1.071 

C2 

-0.274 

-0.342 

-0.176 

-0.239 

-0.102 

-0.161 

C3 

0.021 

0.019 

0.014 

0.012 

0.009 

0.007 

C4 

-0.048 

-0.047 

-0.034 

-0.032 

-0.024 

-0.022 

Cs 

0.010 

0.010 

0.008 

0.008 

0.007 

0.006 

Ce 

-0.061 

-0.058 

-0.039 

-0.037 

-0.026 

-0.023 

Cj / Q^em 

-0.005 

-0.105 

-0.011 

-0.097 

0.035 

-0.081 

Cg/ CXem 

0.086 

0.023 

0.055 

0.014 

0.036 

0.009 

Cg/ CXem 

-1.419 

-0.091 

-1.341 

-0.087 

-1.277 

-0.075 

Cio/ CTem 

0.383 

-0.021 

0.264 

-0.016 

0.176 

-0.011 

C7y 


-0.342 


-0.306 


-0.276 

Csg 


-0.160 


-0.146 


-0.133 
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TABLE II: The values of the Wolfenstein parameters, A, A, p, and rj 


Refs. 

A 

A 

P 

fj 

7 


m 

0.2237 ± 0.0033 
0.819 ± 0.040 ^ 
0.224 ± 0.038 
0.317 ± 0.040 
(54.8±6.2)° 


m 

0.222 ± 0.004 
0.83 ± 0.07 
0.21 ± 0.12 
0.38 ± 0.11 
(62±15)° 


11 

0.2210 ± 0.0020 
0.831 ± 0.022 ^ 
0.173 ± 0.046 
0.357 ± 0.027 
(63.5±7.0)° 


\m 

0.2236 ± 0.0031 “ 
0.824 ± 0.046 * 
0.22 ± 0.10 
0.35 ± 0.05 
(59±13)° 


“Determined from the measurements of |I4d| = 0.9734 ± 0.0008 and |I4s| = 0.2196 ± 0.0026, and the 
error includes scale factor of 1.5 

bA = \Vcb\/X, and \Vcb\ = (40.6±0.8)xl0-3 [||, (41.2±2.0)xl0-3 ||. 

“p = p(l - A72) 

‘^fj = 77(1 — A^/2) 


TABLE III: Values of meson decay constants, form factors, and rj-rj' mixing parameters. 
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TABLE IV: The CP-averaged branching ratios (in the unit of 10“®) of decays Bg PP calculated 
with ^ = mb, and default values of parameters. The data in the column 2~4 are computed with 
A = 0.824, A = 0.2236, p = 0 . 22 , f) = 0.35, 7 = 59°, while the data in the column 5~7 are 
computed with A = 0.82, A = 0 . 22 , p = 0.086, 7 = 0.39, 7 = 78.8°. 


decay 

modes 

NF 

QCDF 

NF 

QCDF 

BR 

BRf 

BRf+^ 

BR 

BRf 

BRf+^ 



8.724 

12.06 

18.81 

7.995 

11.06 

17.25 


0 

0 

T 

0.189 

0.154 

0.200 

0.247 

0.210 

0.277 


0 

T 

0.108 

0.061 

0.071 

0.125 

0.077 

0.091 



0.434 

0.497 

0.717 

0.396 

0.522 

0.777 


0 

0 

T 

— 

— 

0.011 

— 

— 

0.011 



0.052 

0.078 

0.087 

0.058 

0.066 

0.071 


0 

T 

0.055 

0.059 

0.056 

0.061 

0.062 

0.062 



4.570 

7.084 

10.52 

4.087 

6.535 

9.720 



9.190 

10.60 

16.44 

8.401 

9.655 

14.96 


-^p'p' 

4.622 

6.034 

10.65 

4.324 

5.583 

9.846 



9.653 

10.23 

10.44 

7.311 

7.667 

7.700 


^K+K- 

6.949 

9.762 

15.63 

7.807 

10.65 

16.58 


-^7r+7r“ 

— 

— 

0.022 

— 

— 

0.023 
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TABLE V: The CP-averaged branching ratios (in unit of 10“®) of decays Bs PV calculated 
with fi = rrib, and default values of parameters. The data in the column 2~4 are computed with 
A = 0.824, A = 0.2236, p = 0.22, f) = 0.35, 7 = 59°, while the data in the column 5~7 are 
computed with A = 0.82, A = 0.22, p = 0.086, 7 = 0.39, 7 = 78.8°. 


decay 

modes 

NF 

QCDF 

NF 

QCDF 

BR 

BRf 

BRf+‘^ 

BR 

BRf 

BRf+^ 



2.163 

2.812 

4.985 

1.982 

2.581 

4.577 



0.422 

0.630 

2.381 

0.387 

0.577 

2.182 


^K+K*- 

1.818 

2.254 

3.819 

2.598 

3.119 

4.992 


^K-K*+ 

1.262 

1.606 

3.906 

0.863 

1.103 

2.864 


-^ 7 T~^ P~ 

— 

— 

0.002 

— 

— 

0.001 


-^TT~ p^ 

— 

— 

0.002 

— 

— 

0.001 


0 

0 

t 

0.368 

0.367 

0.376 

0.349 

0.350 

0.402 



0.422 

0.428 

0.480 

0.292 

0.314 

0.341 



0.029 

0.057 

0.154 

0.036 

0.070 

0.188 


0 

* 

0 

is 

t 

0.137 

0.068 

0.102 

0.095 

0.052 

0.082 



— 

— 

0.002 

— 

— 

0.001 


^TT^UJ 

— 

— 

0.001 

— 

— 

0.001 


0 

t 

0.086 

0.100 

— 

0.095 

0.098 

— 


0 

* 

t 

0.121 

0.120 

0.230 

0.169 

0.154 

0.300 


0 

* 

t 

0.062 

0.028 

0.038 

0.043 

0.021 

0.028 


0 

t 

0.122 

0.207 

0.223 

0.136 

0.166 

0.174 


0 

t 

0.128 

0.128 

0.125 

0.144 

0.146 

0.148 


^puj 

0.042 

0.043 

0.055 

0.052 

0.031 

0.039 


-^p'oj 

0.044 

0.012 

0.016 

0.055 

0.014 

0.017 


-^pcj) 

0.259 

0.530 

0.456 

0.216 

0.483 

0.417 


^p'cj) 

0.007 

0.266 

0.367 

0.005 

0.238 

0.328 


-^K+p- 

23.69 

24.12 

24.36 

19.02 

19.29 

19.18 


^7r-K*+ 

6.643 

6.949 

6.910 

5.709 

6.008 

6.161 
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TABLE VI: The CP-violating asymmetry parameters a^/ and for Bg PP decays (in the 
unit of percent) calculated with the QCDF approach, with ^ = ruh, and default values of various 
parameters. The data in the column 2~5 are computed with A = 0.824, A = 0.2236, p = 0.22, 
p = 0.35, 7 = 59°, while the data in the column 6~9 are computed with A = 0.82, A = 0.22, 
p = 0.086, 7 = 0.39, 7 = 78.8°. 


modes 

4 

af+“ 


^e+e' 

4 


^e+e' 

f+a 

4+e' 


-0.90 

-0.72 

3.38 

3.43 

-0.98 

-0.79 

3.66 

3.72 

Bg^K^n^ 

-64.72 

-59.92 

-61.67 

-74.32 

-47.54 

-43.46 

-86.01 

-89.93 

Bg^K^sV 

-85.03 

-82.76 

-37.31 

-48.06 

-67.24 

-64.29 

-71.02 

-75.70 

Bg^K^sV' 

54.84 

47.19 

-32.54 

-42.34 

52.27 

43.60 

-39.40 

-47.14 

Bg^TT^TT^ 

— 

0 

— 

-27.71 

— 

0 

— 

-28.25 

Bg^TT^r] 

-16.27 

-15.71 

24.07 

35.66 

-19.20 

-19.39 

26.41 

40.18 

Bg^TT^T]' 

-22.84 

-22.75 

-31.39 

-44.93 

-21.67 

-20.47 

-32.69 

-45.21 

Bg^rjT] 

1.08 

0.86 

1.97 

1.57 

1.18 

0.93 

2.12 

1.69 

Bg^rirj' 

-0.69 

-0.54 

5.14 

5.49 

-0.76 

-0.59 

5.59 

5.98 

Bg^Vj'vi' 

-3.25 

-2.47 

1.22 

1.40 

-3.52 

-2.67 

1.29 

1.50 

Bg^K^K^ 

-6.89 

-5.33 

-40.89 

-33.57 

-6.32 

-5.04 

-40.74 

-33.89 


— 

0 

— 

-27.71 

— 

0 

— 

-28.25 
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TABLE VII: The CP-violating asymmetry parameters a^/ and for Bg PV decays (in the 
unit of percent) calculated with the QCDF approach, with ^ = ruh, and default values of various 
parameters. The data in the column 2~5 are computed with A = 0.824, A = 0.2236, p = 0.22, 


T] = 

0.35, 7 

= 59°, while the data 

in the column 6~9 

are computed with A 

= 0.82, A 

= 0.22, 

P = 

0.086, p 

= 0.39, 7 = 

78.8°. 







modes 

4 

af+“ 


^e+e' 

4 

af+“ 

^e+e' 

/+“ 

Bs- 


-27.38 

-21.76 

68.58 

45.24 

-28.72 

-20.39 

4.69 

-24.52 

Bs- 


38.51 

32.73 

91.14 

87.34 

52.52 

46.18 

79.67 

88.70 

Bs- 


9.29 

6.07 

-76.99 

-75.32 

7.58 

4.99 

-74.30 

-72.78 

Bs- 


— 

0 

— 

99.48 

— 

0 

— 

92.23 

Bs- 


— 

0 

— 

89.01 

— 

0 

— 

39.27 

Bs- 


-20.88 

— 

-13.09 

— 

-21.30 

— 

-14.58 

— 

Bs 

-^pp^ 

-15.06 

-11.51 

38.43 

45.56 

-18.77 

-14.75 

43.57 

52.53 

Bs- 

o 

t 

-25.73 

-30.00 

-52.01 

-60.37 

-22.49 

-25.30 

-51.79 

-59.53 

Bs 

-^puj 

-18.83 

-12.24 

68.45 

74.11 

-26.04 

-17.28 

78.68 

85.81 

Bs 

-^p'oj 

-71.72 

-59.77 

12.78 

38.52 

-63.65 

-57.69 

-47.14 

-27.12 

Bs 

-^p(j) 

6.56 

7.14 

4.68 

4.13 

7.19 

7.81 

5.00 

4.37 

Bs 

-^p'(i) 

10.95 

8.88 

9.91 

9.71 

12.25 

9.93 

10.60 

10.48 
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TABLE VIII: The CP-violating asymmetry parameters Acp (%) for Bg PP decays calculated 
with the QCDF approach, with /r = m;,, and default values of various parameters. The data in 
the column 3~4 are computed with A = 0.824, A = 0.2236, p = 0 .22, p = 0.35, 7 = 59°, while the 
data in the column 5~6 are computed with A = 0.82, A = 0.22, p = 0.086, p = 0.39, 7 = 78.8°. 


modes 

case 

y\cp 

Af+a 

y\cp 

y\cp 

Af+a 

y\cp 


II 

0.17 

0.17 

0.18 

0.18 

Bg^K^TT^ 

II 

-3.24 

-3.86 

-4.41 

-4.59 

Bg^K^sP 

II 

-2.07 

-2.60 

-3.71 

-3.94 

Bg^K^p' 

II 

-1.49 

-1.99 

-1.83 

-2.24 

Bg^TT^TT^ 

II 

— 

-1.38 

— 

-1.41 

Bg^-K^p 

II 

1.16 

1.74 

1.27 

1.96 

Bg^Tl^p' 

II 

-1.62 

-2.30 

-1.68 

-2.31 

Bg^pp 

II 

0.10 

0.08 

0.11 

0.09 

Bg^pp' 

II 

0.25 

0.27 

0.28 

0.30 

Bg^p'p' 

II 

0.053 

0.064 

0.056 

0.068 

Bg^K^TT^ 

I 

-4.42 

-5.09 

-5.91 

-6.91 

Bg^K^K^ 

II 

-2.06 

-1.69 

-2.05 

-1.70 

Bg^TT^TT^ 

II 

— 

-1.38 

— 

-1.41 
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TABLE IX: The CP-violating asymmetry parameters Acp (%) for Bg PV decays calculated 
with the QCDF approach, with /x = m;,, and default values of various parameters. The data in 
the column 3~4 are computed with A = 0.824, A = 0.2236, p = 0 .22, p = 0.35, 7 = 59°, while the 
data in the column 5~6 are computed with A = 0.82, A = 0.22, p = 0.086, p = 0.39, 7 = 78.8°. 


modes 

case 

y\cp 

Af+a 

y\cp 

y^CP 

Af+a 

y\cp 


III 

0.68 

0.27 

0.74 

0.29 


III 

-0.96 

-0.58 

-1.05 

-0.63 

B^^K+K*- 

III 

14.28 

-7.98 

12.02 

-8.50 

B^^K-K*+ 

III 

-13.09 

9.51 

-13.04 

8.71 

B^g^TT+p- 

III 

— 

4.96 

— 

4.60 

B^g^7r-p+ 

III 

— 

4.96 

— 

4.60 

Bg^K^p^ 

II 

3.35 

2.20 

0.16 

-1.27 

Bg^K^uj 

II 

4.64 

4.44 

4.10 

4.54 

Bg^K^gcP 

II 

-3.82 

-3.74 

-3.69 

-3.62 


I 

-63.28 

-51.12 

-81.61 

-63.76 

Bg^TT^P^ 

II 

— 

4.96 

— 

4.60 

Bg^TT^iO 

II 

— 

4.44 

— 

1.96 

Bg^TT^4> 

II 

-0.70 

— 

-0.78 

— 

'^g^pK*^ 

I 

49.72 

30.43 

38.61 

23.34 

Bg^p'K*^ 

I 

-38.13 

-45.57 

-50.43 

-62.38 

Bg^pp^ 

II 

1.88 

2.24 

2.13 

2.58 

Bg^p'p^ 

II 

-2.66 

-3.09 

-2.64 

-3.03 

Bg^puj 

II 

3.37 

3.67 

3.86 

4.24 

Bg^p'uj 

II 

0.46 

1.77 

-2.51 

-1.50 

Bg^pcj) 

II 

0.25 

0.22 

0.27 

0.24 

Bg^p'cp 

II 

0.52 

0.51 

0.56 

0.55 

Bg^K^p^ 

I 

-2.09 

0.79 

-2.61 

1.00 

Bg^TT^K*^ 

I 

0.06 

-5.10 

0.07 

-5.73 
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TABLE X: Penguin-to-tree ratios and bound on 7 , using the default values of various parameters, 


with the QCDF approach. The bound on 7 is from Eq. (^). The data in row a are computed 
with A = 0.824, A = 0.2236, p = 0.22, fj = 0.35, 7 = 59°, while the data in row b are computed 
with A = 0.82, A = 0.22, p = 0.086, f) = 0.39, 7 = 78.8°. 



\Pc/Tc\ 

1 P'K'K/^TVTV 1 

5 ' 

Rkk 

7 

a 

5.359 

0.282 

8.13° 

0.831 

< 61.79° 

b 

5.740 

0.292 

8.13° 

0.961 

< 83.52° 
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FIG. 1: CP-averaged branching ratios of Bg —> K~^K~ (column 1), K^tt^ (column 2), 
(column 3), and (column 4) decays within the QCDF approach versus the angle 7. The 

dashed lines, solid lines, and doted lines correspond to the default values of various theory inputs, 


wherein ms(2GeV) = 90 MeV, 105 MeV, and 120 MeV, respectively. The dot-shades demonstrate 
the uncertainties due to the variations of the GKM matrix parameters A, A, p, fj (row 1), formfactor 
(row 2), {qh, <Ph) (row 3), {qa, (Pa) (row 4), and overall inputs (row 5). 
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FIG. 2: CP-averaged branching ratios of Bg —> 7 r'*“ 7 r“ decays within the QCDF approach versus 
the angle 7 . The legends on the (dashed, solid, and doted) lines are the same as those in Fig. |l|. 
The dot-shades demonstrate the uncertainties due to the variations of the CKM matrix parameters 
A, A, p, fj (a), weak annihilation parameters {qa, (t’A) (b), and overall inputs (c). 



FIG. 3: Penguin-to-tree ratio iPn-Tr/Pn-Trl (a) and Rkk versus 7 . The legends on the (dashed, solid, 
and doted) lines are the same as those in Fig. ||. The dot-shades demonstrate the uncertainties 
due to the variations of overall inputs. 
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